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HIGHLIGHTS 


•  A  simplified  model  is  obtained  from  a  lithium-ion  battery  electrochemical  model. 

•  Fractional  differentiation  is  used  for  a  low  number  of  parameters  in  the  model. 

•  Some  electrochemical  variables  and  parameters  still  appear  in  the  simplified  model. 

•  A  relative  error  less  than  0.5%  on  the  voltage  is  obtained  under  various  conditions. 

•  The  model  simplicity  and  the  accuracy  are  interesting  for  use  in  automobile  BMS. 
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With  hybrid  and  electric  vehicles  development,  automobile  battery  monitoring  systems  (BMS)  have  to 
meet  the  new  requirements.  These  systems  have  to  give  information  on  state  of  health,  state  of  charge, 
available  power.  To  get  this  information,  BMS  often  implement  battery  models.  Accuracy  of  the  infor¬ 
mation  manipulated  by  the  BMS  thus  depends  on  the  model  accuracy.  This  paper  is  within  this  frame¬ 
work  and  addresses  lithium-ion  battery  modeling.  The  proposed  fractional  model  is  based  on 
simplifications  of  an  electrochemical  model  and  on  resolution  of  some  partial  differential  equations  used 
in  its  description.  Such  an  approach  permits  to  get  a  simple  model  in  which  electrochemical  variables 
and  parameters  still  appear. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

More  and  more  drastic  automobile  pollution  norms  and  oil  price 
increase,  lead  car  manufacturers  to  design  new  vehicles  (stop  and 
start,  hybrid,  or  electric).  The  underlying  idea  is  to  reduce  exhaust 
emissions  in  the  built-up  areas,  either  by  stopping  the  internal 
combustion  engine  when  the  vehicle  is  not  moving,  or  by 
substituting  electric  fuel  for  fossil  fuel.  The  latter  may  involve  an 
electric  motor  and  one  or  more  energy  storage  system.  To  keep 
these  new  vehicles  in  good  working  order,  car  manufacturers  must 
integrate  a  reliable  electrical  energy  storage  management  [1],  That 
is  why,  state  of  charge  (SOC)  and  state  of  health  (SOH)  estimators 
must  be  designed.  To  design  these  estimators,  dynamical  models  of 
the  battery  pack  can  be  a  valuable  tool. 
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In  this  paper,  a  dynamical  model  of  a  lithium-ion  battery  is 
proposed.  Many  models  exist  in  the  literature  for  this  kind  of  bat¬ 
teries  such  as  purely  electric  models  [2]  or  fuzzy  models  [3],  An 
extensive  analysis  on  lithium  ion  batteries  modeling  in  proposed  in 
Ref.  [4],  The  originality  here  is  the  way  this  model  is  obtained.  It 
results  from  simplification  of  an  electrochemical  model  that  de¬ 
scribes  the  battery  behavior  using  partial  differential  equations. 
Interest  of  this  approach  is  that  it  is  possible  to  link  the  electro¬ 
chemical  parameters  of  the  battery  to  the  resulting  dynamic  model 
parameters.  This  model  has  thus  an  important  physical  meaning, 
unlike  purely  electric  models  proposed  in  the  literature.  Another 
interest  is  the  introduction  of  fractional  differentiation  that  helps  to 
describe  some  parts  of  the  model  with  a  small  number  of  param¬ 
eters  again  directly  related  to  the  electrochemical  parameters  of  the 
battery.  As  shown  on  several  electrochemical  devices  [5—8],  frac¬ 
tional  differentiation  permits  to  obtain  models  with  a  low  number 
of  parameters.  Such  a  model  can  thus  be  used  for  SOC  or  SOH 
estimator  design. 
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2.  Electrochemical  model 

The  simplified  model  that  will  be  presented  in  the  sequel  is 
based  on  lithium-ion  electrochemical  presented  in  Refs.  [9,10]  that 
results  in  Newman’s  modeling  approach  [11],  This  model  is  a 
pseudo  2D  based  on  a  representation  of  the  cell  such  as  in  Fig.  1. 

In  Fig.  1,  the  electrodes  are  seen  as  an  aggregation  of  spherical 
particles  (2D  representation)  in  which  the  Li+  ions  are  inserted.  The 
first  spatial  dimension  of  this  model,  represented  by  variable  x,  is 
the  horizontal  axis.  The  second  spatial  dimension  is  the  particle 
radius. 

The  cell  is  constituted  of  three  regions  (two  electrodes  and  a 
separator)  that  imply  four  distinct  boundaries  at: 

•  x  =  0:  negative  electrode  current  collector; 

•  xe  ]0;5„[:  region  1,  <5n  thickness  negative  carbon  electrode  (LixC6, 
MCMB---): 

•  x  =  5n  :  negative  electrode/separator  interface; 

•  xe  ]  Sn  ;  5n  +  <5sep[  :  region  2,  5sep  thickness  separator; 

•  x  =  5n  +  <5sep:  separator/positive  electrode  interface; 

•  xe  ]  5„  +  <5sep  ;  i[:  region  3,  5P  thickness  positive  electrode  made 
of  UC0O2; 

•  x  =  L:  positive  electrode  current  collector. 


As  described  in  Refs.  [9,10],  lithium  concentration  cs(x,  r,  t) 
evolution  in  the  spherical  particle  of  radius  r  =  Rs,  is  supposed 
described  by  the  diffusion  law, 


9Cs  _  Ds  9  /  2dcs\  J  ®rlr-0  / 2 

dt  ~  r2  3r  \  dr  J  j  acJ  _  -jLi ' 

1  9r  lr=Rs 

Lithium  concentration  ce(x,  f)  in  the  electrolyte  is  modeled  by: 


9(rece) 

dt 


1  _  fO 

1  r+,-Li 
c  J  • 


with 


dce  I  fice  I 

lx=0  dx  \x=L 


(3) 


(4) 


Charge  conservation  in  the  solid  phase  of  each  electrode  is 
defined  by  the  Ohm’s  law: 


The  two  electrodes  are  assumed  electrochemically  porous.  Re¬ 
gions  1  and  3  are  therefore  constituted  of  a  solid  phase  (electrode 
material)  and  a  liquid  phase  (electrolyte). 

The  cell  is  supposed  supplied  by  a  current  I(t ).  The  cell  voltage 
denoted  Ubatt(f)  is  defined  by  the  relation 

Ubatt  (t)  =  &OM)-<fe(0,t)-|/(f)  (1) 

where  <ps(L,t )  and  <ps(0,t)  are  respectively  the  positive  electrode 
potential  at  abscissa  x  =  L  and  the  negative  electrode  potential  at 
abscissa  x  =  0;  Rf  denotes  the  contact  resistance  and  A  the  elec¬ 
trodes  surface.  Note  that  all  the  parameters  used  in  the  following 
equations  are  defined  in  Refs.  [9,10], 


*e,r(.;> )  -r  =  0, 


(5) 


with  the  following  limit  conditions  at  the  current  collectors 


eff^sl 
-  fix 


eff  d^s  I 

+  fix 


I 

A’ 


(6) 


and  the  null  current  conditions  at  the  separator: 


fix 


90s 

dx 


=  0. 


x=<5n+<5se 


(7) 


If  $e(x,  f)  denotes  the  electrolyte  potential,  charge  conservation 
in  the  electrolyte  is  defined  by: 


Fig.  1.  Pseudo  2D  model  of  a  lithium-ion  cell  9,10], 


l(Keffil^)  +l(K°ff|ln(Ce))  +J'U  =  °’  (8) 

with  the  limit  conditions 


90e  I  90e  I 

dx  L=o  9x  |x=i 


(9) 


The  four  differential  equations  (2),  (3),  (5)  and  (8)  that  describe 
variables  cSie,ce,0s,0e,  variations  are  linked  by  the  Butler— Volmer 
equation 


dsio 


-  exp 


(10) 


In  equation  (10),  jLl  is  induced  by  overvoltage  t),  defined  by  the 
potential  difference  between  the  solid  phase  and  the  electrolyte 
and  equilibrium  thermodynamic  potential  U: 


»?  =  0s  -  0e  -  U.  (11) 

The  equilibrium  potential  U  is  a  function  of  the  solid  phase 
concentration  at  the  spherical  particles.  System  input  is  the  current 
/(f)  in  equation  (6).  System  output  is  the  cell  voltage  given  by 
equation  (1). 


38 


J.  Sabatier  et  al.  /  Journal  of  Power  Sources  262  (2014)  36-43 


3.  From  the  electrochemical  model  to  a  fractional  dynamic 
model 

The  approach  is  here  to  reduce  the  electrochemical  model 
equations  (1)— (11)  to  simpler  equations  taking  into  account 
simplifying  assumptions  and  using  some  approximations.  This  kind 
of  approach  has  been  recently  used  in  Ref.  [12]  that  supposes  small 
signal  current  or  potential  as  excitations  applied  to  the  battery. 
These  excitations  permit  to  simplify  the  electrochemical  model  and 
to  extract  battery  impedance.  Here  small  signal  current  or  potential 
as  excitations  are  not  used.  The  electrochemical  model  previously 
described  has  been  implemented  under  Comsol  Multiphysics  and 
its  parameters  have  been  adjusted  to  fit  a  SAFT  VL7P  battery 
behavior.  This  model  is  then  used  to  speculate  assumptions  and 
approximations  that  are  now  presented,  and  then  to  propose  a 
simplified  model.  The  assumptions  introduced  are  formally  proved 
and  the  approximation  used  will  then  be  validated  through  the 
simplified  model  accuracy. 

3.3.  Assumptions 

•  Assumption  Hi 

Due  to  the  low  resistivity  of  electrode  materials,  it  is  supposed 
that  the  solid  phase  potentials  are  constants  in  the  two  electrodes: 

</>s(X,t)  =  <t>s  _moy(t)-  (12) 

Moreover,  negative  electrode  is  supposed  connected  to  the 
ground.  Negative  electrode  solid  phase  potential  is  thus  null: 

0sn  (t)  =  =  0s_moy(f)  =  0  Vxe[0,<5n],  (13) 

To  illustrate  assumption  Hi  validity,  the  current  profile  of  Fig.  2 
is  applied  to  the  electrochemical  model.  The  resulting  potentials  are 
represented  by  Fig.  3a. 

•  Approximation  Ai 


It  is  supposed  that  electrolyte  potentials  are  constants  in  the  two 
electrodes.  It  is  moreover  supposed  equal  to  its  mean  value  for  a 
given  electrode: 

0e(*> t)  ~  0e_moy(f)-  (14) 

For  current  in  Fig.  2,  Fig.  3b  shows  the  electrolyte  potentials  in 
the  electrodes  and  in  the  separator.  It  permits  to  assess  approxi¬ 
mation  Ai  quality  and  to  check  that  this  potential  remains  close  to 
0.6  V. 


•  Approximation  A2 

Current  density  jLl  is  supposed  constant  in  each  electrode.  It  is 
moreover  supposed  equal  to  its  mean  value  for  a  given  electrode: 

5j 

j'P(t)  =  JPmean(t)/1=  JT  J  jU(t,X)dx  (15) 

'  0 


with  i  =  n,p  (n  =  negative  electrode,  p  =  positive  electrode). 

For  current  in  Fig.  2,  Fig.  4  shows  the  current  density  jLl(f)  in  the 
2  electrodes  and  permits  to  assess  approximation  A2  quality  for  the 
positive  electrode.  Due  to  the  low  impact  of  the  negative  electrode 
on  the  whole  cell  behavior  (this  electrode  only  create  small  voltage 
variations,  aging  been  not  modeled  here),  assumption  A2  is  also 
adopted  for  this  electrode. 

Equation  (5)  integration  for  negative  electrode  leads  to: 


°eff 


Ix=i5„ 


lx=0 


(16) 


or  given  equations  (6)  and  (7) 


(17) 


and  thus  given  approximation  A2  and  for  the  two  electrodes 

JLmeanW  =  ^  =  JpW  «'  =  n>P-  (18) 


•  Assumption  H 2 

Butler— Volmer  equation  defined  by 

JLmean  (f)  =  a^!0i{exp(^Pi(f))  -  exP  }’  (19) 

with 

Vi(t)  =  0sl-(t)  -  Mt)  -  Uj(t)  (20) 

is  supposed  linear  on  the  current  range  defined  by  table  in  Fig.  5  for 
a  SAFT  VL7P  battery. 

This  assumption  is  now  formally  validated  using  the  electro¬ 
chemical  model  reminded  in  section  2.  Even  with  a  high  value 
current  (a  positive  strobe  followed  by  a  negative  one  with  100  A 
each),  rji  variations  remains  close  to  0  (see  Fig.  6). 

Using  the  first  term  of  the  Taylor  expansion  of  the  exponential 
function  in  equation  (19),  a  linear  approximation  of  this  relation  is 
thus: 


Fig.  2.  Current  profile  used  to  validate  assumption  Hi. 
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Negative  electrode 

Separator 

Positive  electrode 

0  0.002  0.004  0.006  0.008  0.01 

Position  x  in  the  electrode  (cm) 


>  2 


0.5 

0 

-0.5 


-  t 

=  0 

t 

=  10s 

-  t 

=  20s 

-  t 

=  30s 

-  t 

=  40s 

-  t 

=  50s 

'  t 

=  60s 

(a) 

Fig.  3.  Solid  phase  (a)  and  electrolyte  (b)  potential  in  all  the  cell. 


(b) 


i(f)  —  asi!0i 


OigF 

RT 


m(t)  - 


acF 

W 


Vi(t ) 


(21) 


with  aa  +  ac=  1.  Using  only  a  graph,  it  is  then  easy  to  check  that  this 
approximation  leads  to  an  error  less  then  0.2%  for  supply  currents 
close  to  250  A  (maximal  peak  current). 

Using  equations  (19)  and  (20),  then: 


-  W) 


iu 

Ji_  mean 


■  «c) 


(22) 


Ee(Ce(x,p)  -  Ce(x,  0))  =  A  (Df  J^(X,p)) 

or  if  Dff  and  re  are  supposed  constant  as  in  Ref.  [10], 

Deff  a2 

Ce(X,p)-Ce(X,0)  =  -j-^Ce(x,p). 

A  solution  of  equation  (24)  is 

Xy/Fe  -Xy/ce 

Ce(x,p)  =  eV°?c ,  +  eV°fc2  +  Ce(x,  0) 


(23) 


(24) 


(25) 


•  Assumption  H3 

Lithium  concentration  in  the  separator  is  supposed  a  linear 
function  of  x  around  the  initial  value  denoted  Ce;(x,0).  As  in  the 
separator  Ju  is  null,  Laplace  transform  of  equation  (3)  is: 


Position  x  in  the  electrode  (cm) 


in  which  Ci  and  C2  are  two  constants  depending  on  limit  condi¬ 
tions.  As  x^/Te/Jof  is  a  small  quantity,  then 
gXyfe/x/D f  —  i  +xv/fj/ and  equation  (25)  becomes 


Ce(x,p) 


V?e 

\fof 


(C1-C2)  +  C1+C2  +  Ce(x,  0) 


(26) 


which  is  the  equation  of  a  straight  line  if  Ce(x,0)  is  constant  or  also  a 
straight  line.  This  thus  validates  assumption  H3. 

This  assumption  validity  is  also  illustrated  by  Fig.  7. 

In  the  separator,  equation  (8)  becomes 


3 

3x 


a 

ax 


<  9 

t(x,  t)  ax 


=  0, 


(27) 


or  after  integration, 

K  ^eiiX,  t)  + 


a 


dx 


cei  (x,  t)  dx 


Cei(X,t)  =  c. 


(28) 


Limit  conditions  (9)  shows  that  constant  C  is  equal  to  0,  and  thus 
equation  (28)  becomes 


K  ax^d(x,t)  +  KceK°  axCd(x,t)  =  ° 


(29) 


Fig.  4.  Current  density  jLl  in  the  electrodes. 


With  Kce  =  1/Ce/(X,0). 
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VL7  P 

VL20  P 

VL30P  1 

Electrical  characteristics  i 

Nominal  voltage  [V] 

3,6  i 

3,6 

3,6 

Average  capacity  1C  after  charge  to  4.0  V/cell  (Ajh] 

7  ! 

20 

30 

Minimum  capacity  1C  after  charge  to  4.0  V/cell  [Ah] 

6,5 

18,5 

28 

Specific  energy  after  charge  to  4.0  V/cell  [Wh/kg] 

67  1 

| 

89 

97 

Energy  density  after  charge  to  4.0  V/cell  (Wh/dpn3) 

131  i 

187 

209 

Specific  power  (1 0s/50%  DOD)  [W/kg) 

1811  ! 

1413 

1136 

Power  density  [10s/50%  DOD)  (W/dm3) 

3526  J 

2974 

2451 

Voltage  limits  ■ 

Charge  (V)  i 

4.JO  (4.1  for  peak) 

Discharge  (V) 

2.5  [2.0  for  cold  cranking) 

Current  limits 

Max  continuous  current  (A)  | 

100  J 

250 

300 

Max  peak  current  during  10s  (A)  i 

250  1 

- -I 

500 

500 

Fig.  5.  VLxP  SAFT  batteries  features. 


Equation  (29)  links  <pei(x,t)  to  ce!(x,t).  Using  linearity  assumption 
H3,  if  one  of  these  variables  is  known  in  the  entire  cell,  and  if  the 
second  one  is  known  on  the  positive  or  negative  electrode,  equa¬ 
tion  (29)  permits  to  deduce  this  second  variable  on  the  remaining 
electrode. 


3.2.  Model  construction 


Analytical  solution  of  equation  (4)  leads  to  the  transfer  function 
that  links  cs j,e  =  cSj(p,r  =  Rs )  to  jp,  given  by: 


^csi  (P±s) 


±:si  ,e(P) 


Cse,i(P) 

J'P(P) 


-Rs 


FD^a,, 


(30) 


Due  to  coth  function,  equation  (30)  cannot  be  easily  used  for 
simulation  or  control  purpose.  A  frequency  approximation  can  be 


used  to  overcome  this  problem.  Bode  diagrams  of  Hcsn,e(p)  in  Fig.  8 
(but  a  similar  one  can  be  obtained  for  Hcspe(p)),  indeed  show  that 
transfer  functions  HCSiie(p)  can  be  approximated  using  the  fractional 
transfer  function 


(31) 


Parameter  1<m  can  be  deduced  from  the  Taylor  expansion  of 
equation  (30).  Gain  J<u  corresponds  to  the  coefficient  of  the  lowest 
degree  term: 


Km  = 


-3 

RsFasi ' 


(32) 


Then  to  impose  a  similar  high  frequency  behavior  to  equation 
(30)  and  to  its  approximation,  a  limit  study  lead  to  define  the  corner 
frequency  wcsei  by: 


Fig.  6.  tin  and  t]p  overvoltage  for  the  current  profile  of  ±100  A. 


Fig.  7.  Concentration  variation  in  the  electrolyte  during  the  first  10  s. 
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Wcsei  =  -rT-  (33) 

Fig.  8  validates  the  approximation  of  transfer  function  Hcsne(p) 
by  Hcsn,e(P)-  A  similar  conclusion  can  be  obtained  with  transfer 
function  Hcsp,e(p) 

Li+  ionic  current  inside  electrolyte  is  modeled  by  equation  (3) 
with  limit  conditions  (4).  Laplace  transform  of  equation  (3)  is 
given  by: 


p,eCe  -  Hm  ce(t,  x)  =  A  (of Ce)  +  LJ±f,  (34) 

Solution  of  this  equation  is  thus: 


ce(p,x)  =  +■ 


(3  -to), -Li  .  Ce(0,x) 


+  C4  exp  - 


y/Pi 


pce 

.0.5(1  p) 


+  c3  exp 


VD~e 


(35) 


where  ce( 0,  x )  denotes  the  initial  concentration.  A  cell  is  made  of  3 
regions  in  which  concentration  variation  is  thus  driven  by  three 
equations  similar  to  equation  (35).  Thus  the  six  constants  C3„,  C3S, 
C3P,  C4„,  C4S  and  C4p  need  to  be  computed  using  the  following  limit 
conditions 


cen(t,x  =  <5  J  =  ces(M  =  S  j 
Cen(t,X  =  d~  +  <5sep)  =  Cep(t,X  =  5~  +  <5sep) 
0Cen(t,X)| 


dx 


=  0 


lx=0 


<  acep(t,X)| 


dx 


=  0 


Den 


9cen(t,  x)  | 

1  ax 


—  Des£p' 


aces(t,X)| 


dX 


lx=dn 


p  aces(t, x) | 


dx 


—  Deod 


lx=5n+5se 


p  0Cep(t,X) 

eptep  ax 


(36) 


x=5„+5se 


The  six  constants  depends  on  Jl'mean(t,x),  and  thus  due  to  re¬ 
lations  (18)  to  the  current  /(t).  Transfer  functions 
tfesn(p,x)  =  Cen (p,x)//(p),  Hces(p,x)  =  Ces(p,x)//(p)  and  Hcep(p,x)  = 
cep(p,x)//(p)  can  thus  be  deduced. 

Bode  diagrams  of  Hcen(p)  in  Fig.  9  (but  similar  ones  are  obtained 
for  Hcep(p))  shows  that  frequency  response  of  this  transfer  function 
is  similar  to  those  of  a  first  order  system.  It  is  thus  possible  to 
propose  the  following  approximation: 


Hcei(P)  =  /H  ■  (37) 

1  (l)ce  i 

System  identification  applied  on  measures  obtained  with  a  50  A 
current  leads  to  the  following  parameters: 

Keen  =  3.7736e~7  mol ■  citT3  and  wcen 
=  io-°-493445  rad  s-1; 


IQep  =  —  1.368e  6  mol  cm  3  and  wCep  =  10  0  862  rad  s  1 . 

Contact  resistance  and  also  initial  solid  electrolyte  interface  can 
be  taken  into  account  by  with  the  resistance  R/in  the  equation: 


Frequency  (rd/s) 

Fig.  8.  Hcsne(p)  and  H^spnpe(p)  Bode  diagram  comparison. 


Ubatt (P)  =  4>sp(L,t)  (38) 

Finally,  negative  and  positive  electrodes  equilibrium  potentials 
are  supposed  defined  by: 


Un(x)  =  8.00229  +  5.0647Z  -  12.578zy2  -  8.6322e-4z-1 
+  2.1765e-5zy2  -  0.46016  exp[15.0(0.06  -  z)\ 

-  0.55364  exp[  -  2.4326(z  -  0.92)]  (39) 


Frequency  (rd/s) 


Fig.  9.  Bode  diagram  of  transfer  function  Hcen(p)  and  its  approximation. 
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Fig.  10.  Simplified  model  obtained. 


Up(x)  =  85.681y6  -  357.70y5  +  613.89y4  -  555.65 y3 

+  281.06y276.648y  —  0.30987  exp^5.657y1150^ 

+  13.1983  (40) 

where  z  and  y  denote  respectively  negative  and  positive  electrode 
solid  phase  stoichiometries. 

Using  all  the  previous  comments,  a  simplified  model  of  a 
lithium-ion  cell  can  be  represented  by  Fig.  10.  A  comparison  of  this 
model  response  with  a  VLxP  battery  response  in  Fig.  11  highlights 
the  model  accuracy.  The  errors  is  always  less  than  0.4%  for  a  current 
varying  from  -100  A  to  100  A.  This  accuracy  also  validates  ap¬ 
proximations  Ai  and  A2 


Fig.  11.  Model  and  battery  voltage  comparison  and  relative  error  at  T  =  -10  °C  and 
Soc  =  50%. 


4.  Conclusion  and  future  work 

In  this  paper  a  model  for  lithium-ion  batteries  has  been  proposed. 
This  is  a  fractional  model  that  stems  from  simplifications  of  an 
electrochemical  model  and  is  the  result  of  this  model  partial  differ¬ 
ential  equations  solution  (solution  of  diffusion  equation  (2)).  Such  an 
approach  permits  to  get  a  simple  model  (low  number  of  parameters 
in  comparison  with  an  approach  based  on  Pade  approximation  that 
lead  to  an  integer  model  with  a  large  number  of  parameters)  in 
which  electrochemical  variables  and  parameters  still  appear  (pa¬ 
rameters  physical  meaning  is  often  lost  with  an  approach  based  on 
model  reduction).  The  obtained  model  permits  to  simulate  the  bat¬ 
tery  voltage  with  a  relative  error  less  than  0.5%.  This  relative  error 
value  has  been  obtained  through  several  validation  tests  [13]  (not 
presented  here  due  to  the  limited  number  of  pages)  under  various 
condition  (state  of  charge,  temperature,  current  values  and  sign). 

The  interests  of  such  modeling  are  multiples.  Since  the  physical 
parameters  of  the  battery  (diffusion  coefficients,  maximum  solid- 
phase  concentration,  ...)  are  contained  in  the  model  parameters 
(see  Fig.  10),  system  identification  can  be  used  to  obtain  numerical 
values  of  these  parameters  from  input-outputs  data.  These  nu¬ 
merical  values  can  be  then  used  to  evaluate  the  battery  state  of 
health.  Also,  some  electrochemical  variables  still  appear  in  the 
obtained  model,  such  as  the  solid  phase  concentration.  Using  a 
state  observer  (or  a  Kalman  filter),  such  a  model  can  thus  be  used  to 
evaluate  battery  state  of  charge.  In  future  papers,  the  authors  will 
develop  these  applications.  Also,  the  authors  will  present  more 
simple  models,  obtained  using  additional  simplifying  assumptions 
and  will  show  how  these  models  has  been  used  to  estimate  state  of 
charge  and  available  power. 

These  works  on  lithium-ion  modeling  and  monitoring  take 
place  in  the  open  lab  “Electronics  and  Systems  for  Automotive” 
combining  IMS  laboratory  and  PSA  Peugeot  Citroen  society. 
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